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PREFACE 


This investigation was prepared during the period 
February - June 1952 at the United States Naval Postgraduate 
School, Monterey, California. 

The subject was suggested to the author by previous 
experience in maintenance of small high-speed diesel 
engines, during which burnt pistons had to be replaced 
on several occasions. 

The author wishes to extend his appreciation to 
Assistant Professor E. E. Drucker for his interest and 


advice throughout the development of the work. 
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SUMMARY 


The objective of this work was to establish by 
annalytical means the temperature distribution in an 
internal combustion engine piston. For this purpose the 
piston was divided into a disk and a barrel, for which 
boundary conditions were expressed for the exterior 
surfaces and a common conical surface. The problem was solved 
in two parts: first, a solution was found for the case where 
the gas temperature was assumed representable by an average 
temperature; this solution was then modified to account 
for a sinusoidal fluctuation of gas temperature. 

The resulting temperature distribution, within the 
limits of the approximations made, is: 


for the disk: 
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CHAPTER I 
INTRODUCTION 


Upon analyzing a piston of an internal combustion 
engine with a viewpoint to understanding its behaviour 
as a heat dissipating member of the engine, it is readily 
realized that the piston absorbs a certain amount of heat 
through the upper surface of its crown. This heat is partly 
stored in the body of the piston until it reaches its steady 
state temperature; at the same time there is heat dissipation 
from the outside surface of the barrel and rings and also 
through its inside surfaces. After the steady state is 
reached, all heat absorbed must be dissipated through these 
surfaces. 

According to 0. L. Adams (1), five percent of the heat 
liberated by the combustion of the fuel mst pass through 
the piston crown surface into the piston. Of this heat, ten 
percent is dissipated by the lower surface of the crown of 
the piston, as stated by J. L. Hepworth (5). 

Radiation from the gases into the piston only takes 
place through a very small portion of the time cycle and 
it amounts, according to B. Pinkel (9), to only ten percent 
of the total heat being absorbed by all the metalic parts. 
forming the container for the hot gases; thus it may be 
seen that the heat absorbed by the piston itself through 
radiation will be a part of this ten percent depending on 
the relative areas of the piston crown wth respect to the 
area of the rest of the combustion chamber. Throughout this 


work it will be assumed that all heat absorption by the 
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piston takes place by convection. Also all heat dissipated 
by the piston will be assumed to do so by a process of 
convection. 

The temperature of the gases throughout the cycle 
varies widely; the temperature curves for an engine may 
be obtained either experimentally by direct measurement or 
annalytically from knowledge of the cyclé under which the 
engine is operating. In this work the temperature of the 
gases will be assumed, for the sake of simplicity, to 
follow a sinusoidal variation about an average temperature. 

The main objective of this work will be to determine 
the temperature distribution throughout the piston. For 
this purpose a set of boundary conditions will be prescribed 
which resemble as closely as possible the actual conditions 


under which the piston operates. 
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CHAPTER ITI 
STATEMENT OF THE PROBLEM 


1. Description of piston. 

A piston with the characteristics shown in figure 1 
was selected for analysis. The crown will be a flat disk 
of uniform thickness 1. The barrel will be a cylindrical 
sleeve of constant thickness b. 

The boundary surface of disk and barrel will be 
surface a-a as shown in figure 1. This surface will have 
properties common to both disk and barrel and this will 
enable some of the unknown quantities entering the questions 
to be found. 

2. Temperature scales. 

Figure 2 shows a sketch of the temperature pattern 
that a particular point of the disk is expected to follow. 
This figure is given mainly to show graphically some of the 
values being used throughout the development. 

Figure 3 serves the same purpose as figure 2 with 
relation to the barrel. 


The temperature of the gases will be assumed to be3 
Tyas = T,, tle 4m 2TH T 


where n will be the frequency of the cycle. 

The temperature Tg surrounding the outside surface of 
the barrel will be assumed constant. 

The temperatures for the disk will be measured from 


the level designated as in figure 2 and will be expressed 
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for any point of the disk as: 
O « t - Tm 


The initial temperature will necessarily be a constant 
and will be designated as Bo. 

The temperatures for points of the barrel will be 
measured in the ) scale of temperatures as shown in figure 


5, where: 


@ +. tet 


The initial temperature $, will be a constant. 
It is realized that the value of Tg will initially be 
equal to t,, but it is impossible with any degree of ease 
to consider into the problem its initial rise from tp, to 
Tg and thus it will be assumed constant and equal to Tg. 
After a solution is reached, the results will be 
reverted to the standard level of measurement of temperatures 
in degrees F. which is designated as the t scale of 


temperatures. 


5. Boundary conditions. 

The following set of boundary conditions has been 
selected: 

(a). At the lower surface of the disk, the rate of 
heat rejection is ten percent of the rate of heat absorption 
at the upper surface. Since the rate of heat conduction is 
a function of the temperature gradient normal to the surface 
considered, and if we assume the coefficient of heat conduction 


to be constant throughout the piston, we may state: 
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There will be a difference in the areas of the upper 
and lower surfaces determined by the relative values of 
Ry and Ry. In this work we will neglect this difference, 
but it is recognized that it could be included by mlti- 
plying the value given as 10 by the ratio of Rq4 to Ry 
squared. 

(ob). At the upper surface of the disk the rate of 
heat absorption from the gases will equal the rate of heat 


conduction from the surface into the disk, or: 
L ta! > ay Cn - 8...) 
J3 re 0 


(c). At the surface a-a the temperature of both the 
Gdisk and the barrel will be equal. Since these temperatures 
will be measured from different levels, this relation may 
not be used directly without introducing undesirable 
constants. 

Therefore it is necessary to realize that at this 
boundary the isotherms from disk and barrel must coincide 
and furthermore they mst be continuous functions of the 
variables involved. This may be stated symbolically: 


along surface a-a; 
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(d). The lower surface of the barrel will be assumed 
insulated, ors 


(e). The inner surface of the barrel will be assumed 
insulated, or: 


a 


Je r=2R, 


(f). At the outer surface of the barrel the rate of 
heat rejection from the barrel into the surrounding medium 
will be equal to the rate of heat conduction from the 


barrel into its outer surface, or: 
ee 3 nz=R 
Jz n=zRo % 


(g). At time = zero, a point: 270, rz=0, will be 
at a temperature 8b, , ors: 


at 
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CHAPTER IIT 
SOLUTION OF PROBLEM 


1. General solution. 

In solving the problem the disk and barrel will be 
treated separately and the solutions matched at the surface 
ana. 

Both the disk and the barrel involve heat transfer in 
a cylinder; the disk may be considered as a flat solid 
cylinder, while the barrel may be considered as a hollow 
cylinder. Under these conditions the best solution of 
Fourier's Law of Conduction of Heat will be that expressed 


in cylindrical coordinates: 


We = 4 yh 4 le 4 lle (1) 
Sr 

Assuming the variables to be separable a solution as 
found in Appendix I iss 


~al(aeb J (2) 


t-:Ee [Camby + Dain by JA. CoreBh Ged] 
From figure 2 it may be seen that the expected curves 


for temperatures will be of the form: 
a7 
G=6, 6° * 
for a particular point on the disk. 
At the same time a sinusoidal variation will be applied, 
which proceeding along the z-axis will be damped and out of 


phase; therefore our solution would have a form such as: 


li 
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°Z accounts for damping. 


where cz is the phase angle and e” 

Examining the previous equation it may be seen that 
the variables z and T may not be separated as assumed in 
the solution found by equation (2), and thus this equation 
will not yield an exact solution. 

A direct solution of the problem was attempted by the 
use of Duhamel's Theorem (3), but it was found impossible 
to match the result obtained for the disk with that for 
the barrel at the surface a-a. 

Under these conditions it was decided to divide the 
problem into one where the gas temperature was assumed to 
be the average temperature T,, and then proceed to super- 


impose on it the effect of the sinusoidal variation in gas 


temperature. 


2 Solution for the disk. 


From equation (2): 
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From boundary condition (b): 


== b bo 


ors k Je M 
mi). 4 = O ..0 


then: 
x Pe ee 


—— 
= 


Gabi or. | kk 


The solution of this equation for b may only be done 
graphically; it will yield a series of values for b which 
will be designated as bye 
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L 
me Sg etter "hes bgt pe hg [EGY mt 
m >I 
From boundary condition (g): 
20 is a) O+ GO. 
Substituting this in equation (5): 
- O, 


13 


nate modaups —_ 
§L 


ve sea, " pes Pty 


adage nomenon 7 
a rir Ata: en r 


3 ae SEONG 


r i «anya pe ik noe inant ly Re RP 
oA 4 an ah y A ; ie Fh , i hy ve ii My “ Ta oe a ‘. ry M i 8 ‘ ¥ 

. “e 4 oy re hie ed Ly a de al +8 HR Petes ti 
' i ; ‘ s a 

i 2 yy ey 

= We 4, a A ba wey! ; ered ed - in 


Wy 
Par i * 
, wes, 4 oe 
mT RP; 
se 
its | 
1 eee ) 7 
t dP Ls 
; pul 
pf gusta 


y t ‘ f bas i ag 
‘a ie rid fs bi > oe , 
ie nd va ne ¥ Mi : ri 


Bh cab us , ied va 2 i ee) be 


i 
4 


wwe : a f 


Therefore: 


sa sable 
axZ. 0, e foot 3+ hanks | Bleu) (6) 


™m ® 


3. Solution for the barrel. 
From Equation (2): 
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Solving this equation graphically for the value of f, 


a set of values fp are found that satisfy this equation, 
and equation (9) becomes: 
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Acceptance of this assumption involves the necessity, 
in order to maintain the equality of the above equation, 


thats: 


where: n= p 
Also, in order that the effect of time on the tempera- 
ture be equal for the same point on the surface a-a as 
approached by either equation (6) or (10), that: 
p= 
Applying this condition to equation (13), we find: 
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This equation is valid for points along the surface ; 
a-a, were: oe « (e, ea 
and should yield the values of Pp. 

Now, the replacement of equation (15) into equation 
(14) will yield an equation in terms of trigonometric and 
Bessel functions of r, valid for values of r ranging from 
Ri to Ry. Since this equation may not be simplified to any 
appreciable extent, it has been considered more advisable 
to replace the values of z and r of a particular point 


along the surface a-a, such that: 
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If now in equation (12) for n = 1 this same point in 
the surface a-a is investigated: 
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Equations (16) and (17) should give the same value 
for Pp. There is no way of proving that this will be so 
except by working out a particular problem. If in so doing 
it is found that the values of Py satisfy simltaneously 
both equations, this will then mean that the assumption made 
that: & = by 
and the derived expressions: 

Tm = Sp 

and: fp=a 
are justified. 

In case that the values of Pp found by means of one 
of the equations mentioned do not satisfy the other 
equation, there is still the possibility of introducing 
this value in the second equation and determining a 
relationship between R4, b, Lo, and 1 such that both 
equations are satisfied and for which the assumption made 
is valid. This would introduce a limitation to the physical 
proportions of the pistons for which this development is 
useful. The relationship just mentioned would take the forms 
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5. Partial solution using Tp. 


On the light of previous conclusions, the solution is 
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6. Superposition of sinusoidal variation. 
In order to superimpose the effect of the sinusoidal 
fluctuation of the gas temperature about a mean temperature 


Tm» the above equation (21) in a simplified form: 
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after a quarter of a cycle, or at f2— we find that: 
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If this were put into equation (21) it would give an 
answer corresponding to the situation where this new 
temperature would have been imposed on the piston since 
time zero. 
To avoid this, instead ofs 
To ain Lint (23) 
use the following expression given by M. Jakob and G. A. 
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where: 
T, is the amplitude of the oscillation at the surface 
of the disk; 


Pole: 


-3 ji accounts for the change in phase along z. 


accounts for damping effects along z, and 


The above expression was given for a thick plate sub- 
jected to conditions similar to those of this problem. Even 
if in this case the disk may perhaps not be considered a 
thick plate, it is felt that this equation will yield a 
fairly good approximation if it is considered that it is 
common knowledge that the effects of fluctuating surface 
temperatures do not penetrate to any great extent beyond 
the surface exposed. In this particular case this effect 
is further minimized by the fact that the fluctuation of 


temperatures is very rapid. 
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From equations (23) and (24) it may be seen that what 
happens at the surface at time zero will have its effect 


at a depth z with a time lag of: 
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The addition of the term expressed in equation (24) 
will affect in equal form both the equations (21) and (22). 
Since for the barrel the values of z will, with the excep- 
tion of its topmost part, be far greater than those for 
the disk, this effect will only be perceptible in its 
upper section. 

Including equation (24) into equations (21) and (22) 
we have as a final result: 
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7. Conclusions. 

It is realized that the results obtained are not rigorous 
but it is felt that they constitute a fairly good approxi- 
mation to the actual temperature distribution in the piston. 

Further work might be done on the subject by solving 
an actual problem in order to check the practicability of 
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APPENDIX I 
GENERAL SOLUTION 


To find a general solution to Fourier's Law of Heat 
Conduction for the case where the temperature should be 
analyzed in the unsteady state for a body best described 


in terms of cylindrical coordinates, we use; 


a 
gle, Lae «aie 
fe Ja 


where: & 


iu 


tay aay TY 
Assume: t= O(7)~ R (2) x Z G3) 


then: 


say: 
= = ao as df eee 
pe: 


Equation (3) may be rearranged as: 


e245 4, oR: 0 


rr R dn 


which is Bessel's equation of the first kind of order zero. 
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Then: R = A Te(ar) +B K% (ar) (6) 


Equation (4) may be rearranged as: 
d*2 
a3? 


which yields: Z- Curby+ D asm by (7 


+b6>2 = 0 


Equation (5) may be rearranged as: 


ae + XQ (at+b*) = 0 


ae (ateb')r 
which yields: 9:-Ee (8) 


But, since: © = OxKx2z 


we find using equations (6), (7), and (8), thats: 
atla* 409F 5 
é<« one + Dainbs [7 Toler) +B Yo(ae)] 
In this equation the constants a, b, A, B, C, D, and 
E must be found so as to satisfy the boundary conditions 


of the specific problem involved. 
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